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Abstract 

In this letter we present a method of constructing dynamical sys- 
tems with any preassigned number of equilibria by adding symmetry 
to another system with at least one equilibrium point. If the resulting 
system is chaotic, we call this procedure chaotification by symmetry 
since the resulting system is chaotic with many symmetrical equilib- 
rium points. 
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1 Introduction 



There are several well-known situations involving computation of equilib- 
ria, such as Nash equilibria and market equilibria. For these two cases, the 
equilibrium is the result of individual agents acting on their own in a non- 
competitive but selfish manner. See [10] for more details. On the other hand, 
the existence of many equilibria in a system makes its dynamics more com- 
plex and shows some special structure. Examples include the well-known 
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multi-scroll attractors [1-2-4-5-6-7 and references therein], chaotic attractors 
with multiple- merged basins of attraction [3], scroll grid attractors [8], and 
2n-wing and n x m-wing Lorenz-like attractors [12-13]. 

In [14] a methodology was presented for constructing some simple chaotic 
systems with any preassigned number of equilibria by adding symmetry to a 
simple 3-D quadratic system with only one stable equilibrium. 

In this letter, we generalize the above method to any m-dimensional sys- 
tem by introducing a singular transformation to generate a system with an 
explicit formula and n equilibrium points. The importance of this result is 
that systems with many equilibria are generally more complex and display 
chaotic attractors with special structure. See for example [9]. 

2 Generating m-dimensional systems with n 
equilibrium points 

In this section, we present our method to generate an m-dimensional system 
with exactly n equilibrium points. Consider the following m-dimensional 

system: 

2-1 = fl %2 j ■■•) x m) 

A = h { X 1, X 2, -,X m ) 



(1) 



X 'm ~ fm ( x li x 2i ■•■> Xm) 

and consider the following coordinate transformation 



U l = (V X 1 + X li) " COS arCCOS 

Uj = Xj, for j — 2, 3, m — 1 (2) 



(V X l + x2 m) 



sm I - arccos ' — — 



where n is the number of desired equilibria. Transformation (2) can add a 
{xjj = 2,3, ...,m- l}-axis rotation symmetry, R {xj: j= 2 ,3,...,m-i}(f n), to the 
original system. 

The first and the last relations of the transformation (2) are deduced from 
the relation xi+ix m = {u\ + iu m ) n in the complex plane, where Ui+iu m is the 
unknown variable. Let X\ + ix m = r (cos 6 + i sin 6) , where r = \/x\ + x^, 
cos 8 = , = = , and sin 9 = —*=£==. Let ui+iu m = s (cos a + i sin a) , where 

X^+X m yj X 1 +X m 
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s = Ju\ + uL, cos a = , u n x , and sin a = , By comparaison, 

we have r = s 11 and a = e+2k7T ; = 0,1,..., n — 1. Here we can choose 



# = arccos I , = J (k — 0). The other cases are similar, but are complex 

in their analysis. Finally, we can obtain the formulas for u\ and u m as in (2). 
On the other hand, we have X\ = s n cos na and x m = s n sinna. Hence we 
have 



Xi 



u ' _ . _ _ _ . 

u\ = x'j = fj (ip l (ui,u m ) ,u 2 ,u 3 , <p 2 (u!,u m )) , for j = 2,3, ...,m- 1 

"' J m Ft,,., J 1 



(y/u\ + u£j cos (^n arccos ^-^t=^ = l Pi (ui,u rn ) 

Xj = Uj, for j — 2, 3, m — 1 (3) 
= (V M i + M m) sin arccos ^ ^ ^ = V?2 (ui,«m) 

By differentiating this equation with respect to £, we have 

( / \du 1 du m du 1 du rn J-' 1 du m \ du 1 ■> m du 1 -' 1 J 

"3 = X 'j = fl (^1 ( U l> U m) ■/" 

gyi 9'-P2 a f2 9<Pi 
dui du m du-i du m 

(4) 

where f k = f k (ip 1 (ui, u m ) , u 2 , u 3 , ip 2 u m )) for e {l,m}. Here we 
use the facts that (ip 1 (ui,u m ))' = + §f^-u' m and Transformation (2) 

is singular and not defined at the points (0, x 2 , x 3 , 0). Hence the two 
systems (1) and (4) are not globally but only locally topologically equivalent. 
Equation (4) is well defined if and only if l^ 1 ^ and S^-f^ 2 - — ^ z ^ x - ^ 0. 

1 v ' J OUl ' OUl OUm OU\ OUm 1 

But these relations are not true for all variables U\ and u m since, for example, 
9ipi ^' Um ^ can vanish for the set of points satisfying tan I n arccos , = 

\u m \ ' 

The equilibrium points of the new system (4) are the real solutions of the 
algebraic equations u'j = for all j = 1, 2, m. Thus we have 

fj (tti, Um) , U 2 , U 3 ,...,<p 2 (til, Um)) = (5) 

for all j = 1,2, ...,m. Assume that the original system (1) has at least one 
equilibrium point (ai, a 2 , a m ). Then let (61, b 2 , b m ) be an equilibrium 
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point of the new system (4). From (3) we have 



a i = + W^j cos ^arccos 

bj = cij, for 



bj = cij, for j — 2, 3, m — 1 (6) 



sm narccos 



Firstly, from the first and the last equation of (6), we have ybf + b^ = 

(V a i + a rnj " ■ Secondly, in order to obtain n values for b±, we can assume 
that cii = 0, (but a m ^ 0, otherwise the transformation is not defined) that 

is cos ( narccos ( — ^-t- I I = T n ( — ) =0, where T„ is the Cheby- 
shev polynomial of the first kind (well defined here since , ^ = £ (—1,1)) 

which has n simple roots (multiplicity here is zero) of the form p k 



b' 



(fc) 

(\ a m\)^ 



cos { ^ 2k 2 n^ ^j ' ^ or all A; = 1, n, in the interval (—1, 1). Hence b\ has n differ- 
ent values of the form b[ k ^ = (\a m \) ™ cos ^ 2fc 2 ~ 1 - >7r j , for all k — 1, n. The last 

equation of (6) and the fact that = arccos I , = I = arcsin 



implies that a m = \a m \ sin n arcsin — . Without loss of generality, 

V \(l a "4)™ J J 

we can assume that a m > 0. Hence b m = a™ sin 7^. Thus the n equilibrium 
points of the new system (4) are of the form 

b[ k) = aicos(^),k = l,...,n 

bj = cij, for j — 2, 3, m — 1 (7) 
6 m = sin ^ 

The above procedure can generate m-dimensional systems with n known equi- 
librium points. From equations (7), we remark that if the original system (1) 
has q equilibrium points, then the new system (4) has nq equilibrium points. 
If the resulting system is chaotic, then we call this procedure chaotification 
by symmetry since the resulting system is chaotic with many symmetrical 
equilibrium points as shown in [14] for some elementary examples. 
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3 Conclusion 



In this letter a method was presented for constructing dynamical systems 
with any preassigned number of equilibria by adding symmetry to another 
m-dimensional system with at least one equilibrium point. The importance of 
this result is that systems with many equilibria are generally more complex 
and display chaotic attractors with special structures as described in the 
current literature. 
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